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The study of nuclear physics using lattice QCD is hindered by an exponentially large signal-to- 
noise problem which is conventionally alleviated by raising the quark masses to unphysically high 
values. We propose a novel form of partial quenching for calculations involving nucleons in which the 
sea quark masses are taken to be smaller than the valence quark masses. It is shown that lowering 
the sea quark masses toward their physical values actually improves signal-to-noise. An optimized 
approach to the physical point in the {nis, m„) plane is proposed, with a full analysis of the cost 
benefit. Improvements in computing time of ~ 10^'-'*"^^ where A is the number of nucleons in the 
system, are shown to be possible. 
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I. INTRODUCTION 

Lattice QCD provides a promising tool for the calcula- 
tion of properties of nuclear systems from first principles. 
In particular, one goal is to use lattice QCD to gain access 
to quantities for which we have little or no experimental 
data, such as the three neutron interaction, necessary as 
input for many nuclear models, or hyperon-baryon in- 
teractions, which may have relevance to the equation of 
state for neutron stars. As improvements in computing 
power and algorithms continue to allow more precision 
in lattice calculations, we are entering an exciting era 
in which the calculation of properties of multiple baryon 
systems is becoming possible, as evidenced by the re- 
cent appearance of the first study of three baryons 1] . 
However, calculations involving baryons with light va- 
lence quarks still suffer from an exponential degradation 
in time of signal-to-noise, resulting in large errors. To 
overcome this problem will require enormous computa- 
tional resources as the number of baryons is increased. 
Creative methods for reducing the signal-to-noise ratio 
(SNR) are necessary if we wish to further explore nu- 
clear physics on the lattice (see, e.g., Q). In this paper, 
we investigate the quark mass dependence of the SNR 
and present a new approach which will greatly reduce the 
computational time associated with these calculations. 

A shared characteristic of most lattice calculations to 
date is the use of unphysically large quark masses. Typi- 
cally this is done because as one lowers the quark masses, 
critical slowing down of the algorithms employed to in- 
vert the Dirac matrix occurs. For calculations involv- 
ing nucleons, the SNR is also improved at larger quark 
masses. Based on the elegant argument by Lepage Q - 
outlined in Sec. II - one can show that the SNR for a cor- 
relator of an operator consisting of interpolating fields for 
A nucleons is approximately ^/j\f e~'^'^^"~2"^^)'^ ^ where 
Mn is the mass of the nucleon, r is the Euclidean time 
separation of source and sink, and Af is the number of 
measurements made. For larger quark masses, explicit 
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chiral symmetry breaking is more severe, and the differ- 
ence Mtv — fm.^ is smaller, thus improving the SNR. 

Raising the quark masses results in systematic errors, 
and extrapolation must be performed to obtain physical 
answers. In an attempt to simultaneously avoid criti- 
cal slowing down and reduce systematic errors, partial 
quenching, in which the valence quark masses are taken 
to be smaller than the sea quark masses, has been em- 
ployed in the mesonic sector. However, it has not been 
clear whether this method would be beneficial in the 
baryonic sector due to a reduced SNR. 

A more careful study of the quark mass dependence 
reveals that an unconventional form of partial quench- 
ing, in which the valence quark masses are taken to be 
larger than the sea quark masses, actually improves the 
SNR. In addition, recent results indicate that it is prop- 
agator production and contractions which consume the 
largest amount of computing time for baryon calculations 
contrary to the mesonic sector. In this paper we 
investigate the quark mass dependence of an array of 
factors affecting the precision of lattice calculations in- 
volving nucleons, and propose that the ideal program for 
approaching the physical limit in baryonic calculations is 
to calculate at physical sea quark mass and extrapolate 
in the valence quark mass only. 



II. SIGNAL TO NOISE ESTIMATES 

Conventionally, hadron properties are computed on the 
lattice by considering correlators of the form G(r) — 
(0|;B(t);B'''(0)|0), where B has some overlap with the 
state of interest. After analytically integrating out the 
fermions, one is left with a new operator, 0(t]A), con- 
sisting of quark propagators from i = to i = r, which is 
a function of the gauge fields. For large time separation, 
the correlator of this object will project out the lightest 
state produced by the operator, G oc e"^""^. From this, 
one can extract the energy of the state (£^o)- 

Here, we outline the arguments presented by Lepage 
to estimate the SNR for correlators calculated on lattices 
with anti-periodic temporal boundary conditions and an 
infinite time extent. Since the correlators are approxi- 
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mated by sampling Af independent gauge configurations, 

1 

GMr) = j^Y.'^{r;A,) (1) 

i—l 

the SNR (TZ) can be computed for large M using the 
central limit theorem, 

a' ^ {\0{r-Ar)-\{0{r;A)}\\ (2) 

At large times the correlator in the first term of will 
project out the lightest state produced by OO^ , so phe- 
nomenological knowledge about the strong interactions 
can be applied to make predictions about the SNR. In 
particular, if B is an operator for producing A nucleons, 
then OO^ will consist of 3 A quark and 3 A antiquark 
propagators from t = to t ~ t. The lightest state pro- 
jected out by this operator will be 3 A pions at rest, so 
we would estimate that our SNR is given by 

7^w VA^e-'^(^^""5™-)-. (3) 

Here, TOtt m^r is the mass of the pion in the partially 
quenched theory where valence quark annihilation is dis- 
allowed. 

Since for physical masses, — ^mj^ ~ 730 MeV, a 
very large number of measurements is required for large 
'''(■'' 1 fni) in order to see a statistically significant 



signal. At shorter time separations, the correlator will 
be contaminated by excited states. In practice, one must 
use a finite time extent, and Eq. [3] has recently been 
shown to give a good approximate upper bound on the 
SNR for T < 2 fm on a lattice with a 4.5 fm time extent 
Above this, backward propagating states must be 
taken into account, and the SNR is expected to be much 
worse. We will concentrate on the range for which the 
Lepage expression holds (Eq. [3]), since it is here that 
measurements are most likely to be made. 



III. SNR IN THE {rus, m„) PLANE 

The main source of the signal-to-noise problem is spon- 
taneous chiral symmetry breaking, which causes the pi- 
ons to be light. Thus, it is expected in general that 
the SNR will improve with heavier quark masses. More 
specifically, it is the ability of OO^ to produce light pi- 
ons which affects the SNR. So one might ask whether it 
is necessary to raise all of the quark masses in order to 
improve the SNR, or only the valence quarks associated 
with the interpolating field. The key to quantifying the 
effect on the SNR of changing sea and quark masses inde- 
pendently lies in partially quenched chiral perturbation 
theory (PQxPT). 

The mass of the nucleon in_SU(2) PQxPT has been 
calculated to order 0{m'^) in j5i]. We have. 
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and contains the ©(m^) terms, presented in [s!]. For the pion we have Q 
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Here, 



m I = Bn(m„ 
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for physical pions, m, 
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), where Bq ^ 2.4 GeV, and 
4 MeV. The parameters 
in the expression for the pion mass are given in Q. Fits 

of the parameters for the nucleon mass up to 0{mq^^) 
were performed to the results presented in [§] for values 
of the pion mass below 400 MeV. The unknown parame- 
ters contained in the 0{rri^) terms for the nucleon mass 



expression were sampled from gaussian distributions and 
used as a measure of the systematic uncertainty. Note 
that, due to the and mwmlg terms, changing the 
sea quark masses affects the nucleon mass at 0{mq), 
while the pion mass is unchanged until O (w,^) . 

Inserting these expressions into Eq. (3), and normal- 
izing with respect to the SNR for physical quark masses. 
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FIG. 1; (Color online) Contour plot in the {mss,mvv) plane 
of the following ratio of SNRs at r = 1.5 fm for = 1 to 4: 
SNR(ms3, m„^)/SNR(mss — m„„ — 140 MeV). Each contour 
represents a factor of two with respect to the dashed curve, 
and darker contours correspond to larger ratios. 



gives the result shown in Fig. 1 for systems with up to 
4 nucleons. We see that the best value for the SNR oc- 
curs at high valence quark mass and low (physical) sea 
quark mass, with improvements over heavy sea quarks of 
about 3 times for a single nucleon, and 100 times for 4 
nucleons. The reason for these improvements can be seen 
by inspecting the mass difference governing the SNR in 
Eq. (3). Keeping the sea quark mass at its physical value 
and raising niw to 400 MeV only raises the nucleon mass 
by about 130 MeV, while the pion mass is raised by 240 
MeV, so that Mjv — ItItt has been reduced from 730 MeV 
in the physical case to 500 MeV. 

Another expression to consider is the time at which 
exponential degradation of signal-to-noise sets in. It was 
shown in [l| that the standard projection of the nucleon 
state onto zero momentum introduces a volume suppres- 
sion to the noise. For small time slices, this suppression 
can dominate, in which case the SNR no longer depends 
exponentially on the number of nucleons. The largest 
time for which this occurs was calculated to be. 
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where the spatial volume should be set by the Compton 
wavelength of the lightest pion in the system, <x . 
For A = 1 we find, using the expressions above for the 
nucleon and pion masses at = 400 MeV, an approx- 
imate 50% increase in tnoise when russ is lowered from 
400 MeV to 140 MeV. 



IV. COST OF CALCULATION 

Based on Fig. 1, in order to optimize the SNR at 
a given valence quark mass, one should lower the sea 
quark masses to their physical values. This suggests that 
an ideal approach to the physical point would be to use 
physical sea quark masses for all measurements, and ex- 
trapolate only in the valence quark mass. This approach 
has the added benefits of requiring a single set of gauge 
field configurations to be produced for all measurements, 
as well as reduced systematic uncertainties. 

To determine whether this approach will be beneficial 
in practice, one must include the cost of gauge field and 
propagator production. There are many factors involv- 
ing the quark masses which affect computation time, in- 
cluding the cost of inverting the Dirac matrix, volume 
requirements, and number of independent sources per- 
mitted per gauge configuration. 

The noise, including volume effects, can be approxi- 
mated by [l| 
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where Lj is the temporal extent, Zi is the overlap onto 
the ith state, and Afsrc = J^oi^MNY L^Lt is the num- 
ber of independent measurements which can be made on 
a single configuration. Based on the results in we 
have chosen the normalization Ao such that at niss = 
rriyy = 390 MeV we have 200 sources per configuration. 
To determine the cost of achieving a fixed SNR at a given 
quark mass, we found the number of gauge configurations 
necessary, Afg, then multiplied by the cost of producing a 
single gauge configuration plus the measurements made 
on that configuration. To explore the effects of different 
cost functions we chose both domain wall and staggered 
fermion actions, and domain wall fermion propagators. 
The cost functions used are given in the Appendix. 

Figure 2 shows the ratios of the total costs to perform 
calculations for A = 1 to 4 nucleons for a fixed 10% sta- 
tistical error at niyy — 400 MeV for two different values of 
sea quark mass: niss — 400 MeV versus niss = 140 MeV 
(circles) . We also compared two separate approaches for 
extrapolation to the physical point to determine whether 
our proposal will be beneficial as one lowers the valence 
quark masses. First, we calculated the cost of producing 
gauge field configurations and propagators for equal sea 
and valence quark masses at 8 different values between 
400 and 250 MeV, again at a fixed 10% statistical er- 
ror. Then, we compared this with the cost of producing 
a single set of gauge configurations at physical sea quark 
mass, as well as propagators on these configurations at 
the same 8 values of valence quark masses between 400 
and 250 MeV (squares). 

For this plot we chose a — 0.093 fm for the lattice 
spacing and domain wall fermion cost functions for both 
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FIG. 2: (Color online) Log plot of the following ratio of costs 
at fixed SNR for calculations at different values of quark 
masses with domain wall fermions: Cost(mi,s = m„i, = 400 
MeV)/Cost(mss = 140 MeV, m„„ = 400 MeV), as a function 
of the number of nucleons (circles). Also plotted is the ratio of 
costs for eight measurements using two different approaches 
to the physical point: Cost(250 MeV < niss = < 400 
MeV) /Cost (m,, = 140 MeV, 250 MeV < m„„ < 400 MeV) 
(squares) . 



culations in which theoretical tools for extrapolation to 
the physical point are less well developed. In addition, 
extrapolating to the physical point in only one quark 
mass will greatly reduce the proliferation of fit parame- 
ters usually associated with partial quenching. However, 
there may still be extra non-analytic terms introduced, 
which in principle will require more measurements to ac- 
curately determine their coefhcients. 

A further consideration not addressed in this work is 
the possibility that lowering the sea quark mass will de- 
crease the gap between the ground state and the first ex- 
cited state. This could force one to make measurements 
at larger time separations in order to extract the ground 
state signal. This issue is highly dependent on the system 
one wants to consider, and can be improved by optimiz- 
ing the interpolating field used, as well as the fitting tech- 
nique. However, as observed in eliminating excited 
state effects from calculations involving multiple baryons 
can be difficult even for large sea quark masses. Because 
present day techniques may not be sufficient to extract 
many nuclear observables of interest, further study in 
these areas, particularly as calculations continue to move 
closer to the physical point, will be necessary. 



sea and valence quarks. For A > 1, we found negligible 
sensitivity to changing either the lattice spacing or the 
cost function. Some sensitivity to the cost function was 
found for A — 1. All error bars reflect the uncertainties in 
the HBxPT expressions for the masses, except for A ~ 1 
where the uncertainty due to the choice of cost function 
is also added in quadrature. 

V. DISCUSSION 

For all but the single baryon case, it is clearly benefi- 
cial to calculate nucleon properties at physical sea quark 
masses. From the linear dependence of the log plots in 
Fig. 2 we see that improvements of ~ 10^(^-1) for a 
single calculation and ^ 10'^"^ for a full approach to the 
physical point can be expected. It is still unclear whether 
our method would be beneficial for single nucleon calcu- 
lations unless the gauge field configurations were already 
available. 

Note that we are only considering here a fixed statis- 
tical error. Lowering the sea quark masses will also help 
reduce the systematic errors associated with unphysical 
quark masses. This will be particularly beneficial for cal- 



VI. CONCLUSIONS 

We have shown that to optimize the SNR for nu- 
cleons in terms of valence and sea quark masses, one 
should choose heavy valence quark masses and physical 
sea quark masses. This choice not only improves the SNR 
as compared to the standard choice of heavy valence and 
sea quarks, but should also produce results which are 
closer to the physical case of interest, thus reducing sys- 
tematic errors. These improvements become more sig- 
nificant as the number of baryons is increased, possibly 
making previously intractable calculations realistic in the 
near future. 
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TABLE I: List of parameters for the given cost functions. 



A 


7.98 X 10"^ 


B 
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Ci 


0.01021 


Ca 
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Cm 


1 


cl 


4 


Ca 


4 
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Appendix: Cost Functions 

Example cost functions for domain wall fermion prop- 
agators, domain wall fermion gauge field configurations, 



and staggered fermion gauge field configurations were 
taken from [§], and respectively. Tlie parame- 
ters used in this work are given in Table 1. 



Domain wall propagators: 



Cost cx [A + ) L^LtL^ 

Ttljyy 



Domain wall gauge configurations 

Cost oc 

im / V fm 



(A.l) 



/MeV\ / fmy 
V "iss / V a / V J 



Staggered gauge configurations: 

/20 MeYY"- f Lx Y'' /O.l fm' 
Cost cx ' ' ' ' ' 



\ rus ) \3 fm 



(A.3) 



Here, is the sea quark mass, and rriK is the mass of 
a kaon containing a light sea quark. 



